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Abstract 

We solve the equation of the equilibrium of the gravitating body, with 
a polytropic equation of state of the matter P = Kp 1 , with 7 = 1 + 1/n, in 
the frame of the Newtonian gravity, with non-zero cosmological constant 
A. We consider the cases with n — 1, 1.5, 3 and construct series of 
solutions with a fixed value of A. For each value of n, the non-dimensional 
equation of the static equilibrium has a family of solutions, instead of the 
unique solution of the Lane-Emden equation at A = 0. The equilibrium 
state exists only for central densities po larger than the critical value p c . 
There are no static solutions at pa < p c - We investigate the stability of 
equilibrium bodies in presence of nonzero A, and show that dark energy 
decease the dynamic stability of the configuration. We apply our results 
for analyzing the properties of equilibrium states for cluster of galaxies in 
the present universe with non-zero A. 



1 Introduction 

Detailed analysis of the observations of distant SN la (Riess et al., 1998; Perl- 
mutter et al., 1999) and of the spectrum of fluctuations of the cosmic microwave 
background radiation (CMB), see e.g. Spergel et al. (2003), Tegmark et al. 
(2004), have lead to conclusion that the term, representing "dark energy" con- 
tains about 70% of the average energy density in the present universe and its 
properties are very close (identical) with the properties of the Einstein cosmolog- 
ical A term. This value of the cosmological constant was anticipated by Kofman 
& Starobinsky (1985), basing on the analysis of the existing upper limits for the 
microwave background anisotropy, see also Eke et al. (1996). In papers of 
Chernin (see reviews 2001 and 2008) the question was raised about a possible 
influence of the existence of the cosmological constant on the properties of the 
Hubble flow in the local galaxy cluster at close vicinity of our Galaxy. Basing 
on the observations of Karachentsev et al. (2006), it was concluded by Chernin 
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(2001, 2008) that the presence of the the dark energy (DE) is responsible for 
the formation of this Hubble flow. This conclusion was criticized by Lukash & 
Rubakov (2008). 

The importance of the DE for the structure of the local galaxy cluster (LC) 
depends on the level of the influence of DE on the dynamic properties. In 
particular, it is necessary to check, if the LC may exist in the equilibrium state, 
at present values of DE density, and the LC densities of matter, consisting of 
the baryonic and dark matter (BM and DM). 

Here, we construct, in Newtonian approximation, exact equilibrium configu- 
rations of the sclfgravitating objects in presence of DE. We consider a polytropic 
equation of state, having in mind, that very extended configurations, like LC, 
where the dark matter have not reached the state of a full thermodynamic equi- 
librium, and effective temperature could be not a constant over the cluster. In 
presence of DE (A) the polytropic configurations are represented by a family 
of configurations, instead of the single model for each polytropic index n. The 
additional parameter (3 represents the ratio of the density of DE to the matter 
central density of the configuration. For values of n = 1, 1.5, 3, corresponding 
to the polytropic powers 7 = 2, 5/3, 4/3, we find the limiting values of (3 C , 
so that at /? > /3 C there are no equilibrium configurations and only expanding 
cluster, may be with a Hubble flow, may exist. We apply these results to the 
LC using the observational data and different suggestions, connected with DM 
distribution. Note that the structure of the isothermal configurations in a box, 
in presence of DE was investigated by de Vega & Siebert (2005): in this paper 
the authors were interested in the problem of thermodynamic stability of such 
configurations, following the works of Antonov (1962) and Lynden-Bell & Wood 
(1968). We have derived a virial theorem for the equilibrium configurations in 
presence of DE, and investigate an influence of the DE on the dynamic stabil- 
ity of the equilibrium configurations, using an approximate energetic method. 
It is shown that DE produce a destabilizing effect, contrary to the stabilizing 
influence of the cold dark matter (McLaughlin & Fuller, 1996). 



2 Main equations 

Let us consider spherically symmetric equilibrium configuration in Newtonian 
gravity, in presence of DE, represented by the cosmological constant A. In this 
case, the gravitational force F g which a unit mass undergoes in a spherically 
symmetric body is given by the formula (Chandrasckhar, 1939) 



F * = —3- + -> W 



Gm Ar 

~ + T 

where m is the mass of the matter inside the radius r, which is connected with 
the matter density p by the continuity equation 



dm o 

— =47Tflr 2 . (2 
dr 
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The density of DE p v is connected with A as 

Then, the equilibrium equation of the selfgravitating body is given by 

IdP Gm Ar 
p ar r A 3 

Now, let us consider a polytropic equation of state 

P = Kp J , with 7=1 + -. (5) 

n 

By introducing the nondimensional polytropic variables £ and 8 n so that (Chan- 
drasekhar, 1939) 

r = at; and p = p 6>^, (6) 

we obtain, from (J2]) and (j4]), the equation describing the structure of the poly- 
tropic configuration in presence of the DE as 

Here P is the pressure and po is the matter central density, a is the characteristic 
radius given by 

2 (n + l)K i_i 
« =-i^-^ ' (8) 
while /? = A/iirGpo — 2p v /p represents the ratio of the DE density to the 
central density of the configuration. At j3 = we obtain from the well 
known Lane-Emden equation (Chandrasekhar, 1939). 

Note, that the generalization of Lane-Emden equation in presence of A was 
first obtained by Balaguera-Antolmez et al. (2007). 



3 The virial theorem 

Let us find the relations between gravitational e g and thermal £th energies, and 
the energy £a, of the interaction of the body with DE. The gravitational energy 
of the selfgravitation is written as (Landau & Lifshitz, 1980) 



_ f M mdm t f r , 2) , r 2 (IdP Ar\ 

(9) 

Here m is written using the equilibrium equation (|4]). For the adiabatic stars 
with a polytropic equation of state we have the relations 
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P n + 1 

P E = nP 1 I = E+ — = -^—E, (10) 
p n 

where E and I are thermal energy and enthalpy per mass unit. After some 
transformations, similar to ones made by Landau & Lifshitz (1980), we obtain 
the following relations between the energies 

3 n In 

£g = Eth + 2SA, Eth — -7T £ g + IT £ A, (11) 

Th O O 



£tot = £th + £ g + £\ = — 5 — £ g H 5 — £a = — - — £th + 3ea- (12) 



3 — n 2n + 3 n — 3 
— £ 9 + -g-e A = — 

where £a = — 47r ^ p " f Q M r 2 c?m = — — J Q M r 2 dm is defined below in (|4"0"]l and the 
additive constant in the energy definition of ea is chosen so that ea = at 
A = and M = 0, being e t /j = J Q dm. Now, let us make an additional 
transformation of the expression for s g . Then, by partial integration, we obtain 

f M mdm GM 2 G f R m 2 

We can transform the last integral by using ([9]), (| 10[) and making partial inte- 
grations. We have 

G f R m 2 1 f R 2 (\ dP Ar\ m , 1 f R mdP , A f R 

-dr = / r — -dr = / — dr-\ — / mr dr 







: 2 Jo \p dr 3 / r 2 2 J p dr 6 J 

L f M ldm+ —MR 2 - — [ M r 2 dm= — e th + A M i? 2 + ^. (14) 



2 7 12 12 7 2n 12 2 

Then, by using (fTTj) , (|14l) . we obtain from (fTB"]) the following relations 

3 GM 2 A - 2n + 5 



Mi? 2 - ea, (15) 

5 — n 

Eth= + , Affl 2 + -^£ A . (16) 

' ' 5 — n 



5 


- n i? 


2(5 - n) 




n GM 2 


nA 


■5 


-n R + 


6(5 -n) J 


16 


| in (HU), we 


get 


n 


- 3 GM 2 


(n - 3)A 


5 


-n R + 


6(5 - n) 



Etot = t — + — MR + EA- (17) 

We can calculate E tot for some particular cases. For n = 3, 1, and we have, 
respectively. 



1 GM 2 1 4 _ 2 1 3 GAf 2 1 ( _2 
£tot = Se a , Etot = --——-—AMR z + -e a , £tot = -7— ^ TK KMR ■> 

2 R 12 2 5 it 10 



(18) 



4 



The Lane-Emden model with n = 5, has an analytical solution with finite 
mass M, finite values of the energies and infinite radius R, so that must be 
(5 — n)R —¥ const, if n —¥ 5. In presence of DE the finiteness of values of all 
types of energies require that 

(5 — n)R -> const and ea -> -^-MR 2 , if n -> 5. (19) 

The Lane-Emden solution (without DE) at n — 3 has the zero total energy, 
with a given finite radius and corresponds to a neutral equilibrium. So, the 
knowledge of the total energy of the configuration permits to find the boundary 
between dynamically stable (n < 3, Etot < 0) and unstable (n > 3, Etot > 0) 
configurations. In our case the virial theorem does not permit to do it, because 
the value of ea is not properly defined, while the presence of DE in the whole 
space does not permit to chose, in a simple way, a universal additive constant of 
the energy. Therefore, in spite of the result from (fig)) , where e to t = 3ea < at 
n — 3, in accordance to the stability analysis made in section 4.3, the polytropic 
solution at n — 3 in presence of DE becomes unstable. 

The virial theorem in presence of A was investigated by Balaguera-Antolincz 
et al. (2007). The final result in this paper was presented in the form of integrals, 
instead of simple algebraic relations obtained above. 



4 Equilibrium solutions 

4.1 The case n=l 

At n = 1 the equation (J7J) is linear and has an analytic solution. Following 
Chandrasekhar (1939), let us introduce a new variable x = After the 

substitution the equation (JT]) is transformed into 



, -X + £3. (20) 



fx 

At j3 = we have a well known analytic polytropic solution (Chandrasekhar, 
1939) 

* = ^, (2D 

satisfying the boundary conditions in the center #i(0) = 1, 6*^(0) = 0. At 
nonzero j3, the solution satisfying the boundary conditions in the center is writ- 
ten as 

0i = (1-j8)^+j8. (22) 

The radius of the configuration is defined by the first zero of the solution (l2~Tj) , 
or (f2"2"| . For the Lane-Emden solution ([2"TT) the outer radius corresponds to 
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£,out = 7T and, for the solution of (|2"2"| with DE, the radius of the configuration 
is determined by the transcendental equation 

(1-/3)^^+^ = 0. (23) 

This equation has real solutions only at (3 < (3 C , so that at the outer boundary 
not only 0\ = 0, but also 9[ = for /3 = j3 c . It follows from 



(l-fll^-^)- (24) 



cos£ sin£ s 

therefore the parameters /3 C and £ ut,c of the limiting equilibrium solution in 
presence of DE are determined by the algebraic equations 



(1 - ftc) ° u '° + C = 0, tan Co«t,c = fout.c, < Co«t,c < — • (25) 

The density distribution for the equilibrium configurations with values (3 = 
0, (3 = 0.5/3 c , (3 — f3 c is shown in Fig.l. The nonphysical solution at j3 = 1.5/3 C , 
which has not an outer boundary, is given by the dash-dot line. The nonphysical 
parts of the solutions at f3 < C , behind the outer boundary, are given by 
dash lines. At large £, these solutions asymptotically approach the horizontal 
line Q\ — (3. Numerical solutions of equations (|25j) and (|23|) give j3 c = 0.178 
and £ oui = 7T, 3.490, 4.493, for = 0, /3 = 0.5/3 c = 0.089, /3 = (3 C = 0.178, 
respectively. 



4.2 Some general relations 

The equilibrium mass M n for a generic polytropic configuration which is solution 
of the Lane-Emden equation is written as 

f R 

M n = 4tt pr 2 dr = 4tt 
Jo 

Using equation J?]), the integral in the right site may be calculated by partial 
integration, what gives the following relation for the mass of the configuration 

Reminding (|8]l. we can write the expression for the mass in the form 

Let us note, that here 6 n (£) is not a unique function, but it depends on the 
parameter (3, according to ([7]). For the limiting configuration, with (3 = f3 c , we 
have on the outer boundary £ ou t 



_ f2 | ^ 
^out \ J/- 



(16, 



(28) 



G 



o, 



d9„ 



16™ 



0. 



(29) 



Therefore the mass M n ,Um of the limiting configuration, with account of ([6]), is 
written as 



A-KpQcOi 



— -^-fout ftcPOc — -^-r out " p c , 



(30) 



so that the limiting value /3 C is exactly equal to the ratio of the average matter 
density p c of the limiting configuration, to its central density po c 



POc 



(31) 



For the Lane-Emden solution (with /3 = 0) we have po/p — 3.290, 5.99, 54.18 
for n = 1, 1.5, 3, respectively. Let us consider the curve M(po) for a constant 
DE density p v — A/87rG. For plotting this curve in the nondimensional form, 
we introduce the characteristic density p c h and write the expression for the mass 
in the form 



M n = 4tt 



(n + l)K 



4ttG 



3/2 



Pch 



_p2 I dOn 



Bf 3 



(32) 



where po = po/p c h is the nondimensional central density of the configuration. 
We introduce also the nondimensional mass of the confuguration M n as 



M n = 4tt 



(n + l)K 



-i 3/2 



4ttG 



Pll *M n , Mn=P% 



-f 2 

^out 



+ 



(33) 

Now we can plot the nondimensional curve M n (po), at constant p v = f3po/2 
(from the definition of B). It is convenient to construct the curve starting from 
the model with po = 1 and f3 <§; 1, very close to the Lane-Emden solution, and 
than following the sequence by decreasing the central density po with taking 
(3 oc l/po (being constant p v ) until the value (3 = fi c . For n = 1, using ([24]) . we 
have 



Mi = /5q 



(1 - P){s\nt] out - t] out cos £ out ) + 



P£,out 



PoP = Pi, 



const. 



(34) 

The behavior of Mi(/5o)|a (at constant p v ) is given in Fig. 2 for different values 
of B in = 0, /3 m = 0.5/3 c , /3 m = P c , for which Mi = n, 3.941, 5.397 at p = 1. 
We note that for /3i n = j3 c there are equilibrium models only for po > 1. 
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4.3 The case n=3. Dynamic stability 



The equilibrium equation for this case is 



1 A 



+ p. 



(35) 



The mass of the configuration is given by 



K 



-i 3/2 



Ma, M 3 = -e out (^ 



) 



pa 



M 3 = 4tt 



ttG 



+ 



3 



(36) 



out 



The Lane-Emden model (/? = 0) has a unique value of the mass, independent 
on the density (equilibrium configuration with neutral dynamical stability) . At 
(3^0 the dependence on the density appears because the function 9 3 is different 
for different values of p and, along the curve M 3 {po)\^, the value of /3 is inversely 
proportional to po (like in the case n = 1). 

The density distribution for equilibrium configurations with ft = 0, P = 
0.5/3 c , P = f3 c is shown in Fig. 3. The nonphysical solution at P = 1.5/3 C , which 
has not an outer boundary, is given by the dash-dot line. The nonphysical 
parts of the solutions at /? < f3 c , behind the outer boundary, are given by the 
dash lines. At large £, these solutions asymptotically approach the horizontal 
line 9 3 — /3 1 / 3 , with damping oscillations around this value. The numerical 
solution of the equation ((SSJ) gives (3 C = 0.006, £, out = 6.897, 7.489, 9.889, for 
P = 0, p = 0.5/3 c = 0.003, P = p c = 0.006, respectively. In Fig.4 we show the 
behavior of M 3 (po)|a, I0r different values of /3j„ = 0, Pi n = 0.5/3 c , Pi n = P c , for 
which M 3 = 2.018, 2.060, 2.109, at p = 1, respectively. 

The behavior of M3(po)\a in Fig.4, showing a decreasing mass with the 
increase of the central density, corresponds, for an adiabatic index equal to 
the polytropic one, to dynamically unstable configurations, according to the 
static criterium of stability (Zcl'dovich, 1963). When the vacuum influence is 
small, it is possible to investigate the stability of the adiabatic configuration 
by the approximate energetic method (Zeldovich & Novikov, 1966; Bisnovatyi- 
Kogan, 2001). For n — 3, at po ^> p v , the density in the configuration is 
distributed according to the Lane-Emden solution p = Pq0 3 (^). In this case 
we may investigate the stability to homologous perturbations, changing only 
the central density at fixed density distribution, given by the function #3. Let 
us first calculate the newtonian gravitational energy of the configuration at 
presence of the cosmological constant (DE). For spherical configurations, the 
Poisson equation for the gravitational potential , in presence of DE is given 



The gravitational energy of a spherical body e g , connected with its selfgrav- 
ity, is written as (Landau & Lifshitz, 1980) 



by 



A(p : 
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r M m dm ( r 

e g = -G , m = 4ir pr 2 dr, M = m{R), (38) 

Jo r Jo 

where R is the radius of body. Here, differently from the gravitational potential 
(p where we have the normalization ip = at r = oo, for (p\ with uniform density 
p v this normalization is not possible. Then we can choose tpA — at r = as 
the most convenient normalization. This choice, together using the equation 
(|37|) . leads to the following expression 

knGpv 2 

¥>a = g — r . (39) 

Consequently, the gravitational energy e\, connected with the interaction of 
the matter with DE, is given by 



E\ = / ip\dm — — / r 2 dm. (40) 



o 



In order to analyze the effects due to the presence of DE on the dynamical 
stability of spherical systems with a polytropic equation of state, it is necessary 
to consider configurations close to the polytropic (adiabatic) equilibrium solu- 
tion at n = 3 (and (3 = 0), where the turning point of stability is expected. In 
this case, the presence of DE docs not affect significantly the gravitational equi- 
librium because the nondimensional density profile of the configuration keeps 
constant during the adiabatic perturbation (homologous displacement), so that 
we can neglect the term depending on A in the Lane-Emden equation and the 
nondimensional quantities £ ou t and #3(£) of the unperturbed polytropic solution 
at n = 3 can be used in calculating the gravitational energy e g connected to 
the selgravity. Therefore, also the terms depending on A which appear in the 
expression (fT5|) of e g can be deleted and the gravitational energy e* becomes 



^ f M mdm 4ttG Pv f M 2j 3 GM 2 A f M 2j 

e*=e g +e A = -Gj o — —J rdm = --—--J r dm, 

(41) 

where, taking into account the nondimensional variables ((6]), the energy e\ can 
be written as 

£A = --j^ r 2 dm = --^po« 5 A eU 4 dC (42) 

Here, the last integral in (|4"2"]l has been calculated by Bisnovatyi-Kogan (2001) 
so that 

tout 

3 3 £ 4 d£ = 10.85. (43) 

o 

In the analysis of the dynamical stability, let us consider the total energy e 
of the configuration. Then, in addition to the Newtonian gravitational energy 
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e g and the energy e\, we must also take into account the contribution of the 
thermal energy E t h, corresponding to a specific energy density E (energy per 
unit mass), and we may include a small correction eqr due to general relativity 
(Bisnovatyi-Kogan, 2001). We obtain 

£ = Eth + £g + £A + £GR = (44) 

Edm - 0.639GM 5/3 pl /3 - 0.104AM 5 / 3 p~ 2/3 - 0.918 5 p 2 /3 , 



o 

where we used the following relations for the polytropic configuration with n = 
3: 

M l/3 -1/3 

6^ = 6.897 and R = aUt = . ™ a ■ (45) 

The equilibrium configuration is determined by the zero of the first derivative 
of £ over po, at constant entropy S and mass M, while the stability of the 
configuration is analyzed through the sign of the second derivative: if positive, 
the configuration is dynamically stable, if negative, the configuration is unstable. 

1/3 

It is more convenient to take derivatives over p than over po- Then, we 
have 



d£ = 3p - 4 / 3 f 1 P J m -0.639GAf 5 / 3 +0.208AA/ 5 / 3 Po 1 -1.84^ff p y 3 = 
Jo 9{m/M) ' c z 

(46) 



for the equilibrium configuration and 



32e ^ /3 f H)^-°-- 5 'v«-- G2M7/3 



/3 



2 



(47) 

for the analysis of the dynamical stability. Here 



are the adiabatic index 7 at constant entropy S* and the nondimensional function 
(j), remaining constant during homologous perturbations, respectively. It follows 
from (14"?| that DE input in the stability of the configuration is negative and sim- 
ilar to the influence of the general relativistic correction (Chandrasekhar, 1964; 
Merafina & Ruffini, 1989). On the contrary, the presence of the dark matter, as 
a background, has a positive influence, increasing the dynamic stability of the 
configuration (McLaughlin & Fuller, 1996; Bisnovatyi-Kogan, 1998). Therefore, 
an adiabatic star with a polytropic index 4/3 which is in the state of a neutral 
stability in the pure Newtonian gravity, becomes unstable in presence of DE. 
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The presence of DE decrease the the level of relativity of a warm dark matter 
(WDM), at which it is still possible to form gravitationally bound configurations 
from WDM. 

Note that dynamic stability of pure polytropic models was investigated by 
Balaguera-Antolinez et al. (2007), using static criterium of stability based on 
M(p c ) dependence (f3"2"|) . Our criterium (|47p is valid for any equation of state 
P(p), for configurations near the point of the loss of stability. Influence of DE 
on the position of the bifurcation point for ellipsoid-spheroid transition was 
investigated by Balaguera-Antolinez et al. (2006). 



4.4 The case n=1.5 

The equilibrium equation for this case is 



The mass of the configuration is written as 



M 3/2 = 4tt 



5K 



inG 



3/2 



1/2 



Pch M 3/2 



M 3/2 = 



.1/2 

Po 



2 / M 3/2 

tout y d £ 



out 



(50) 

The density distribution for equilibrium configurations with ft — 0, ft — 0.5ft c , ft — 
ft c is shown in Fig. 5. The nonphysical solution at ft = 1.5ft c , which has not 
an outer boundary, is given by the dash-dot line. At large £, these solu- 
tions asymptotically approach the horizontal line 8 3 / 2 = /3 2 / 3 . The numeri- 
cal solution of equation (09) gives ft c = 0.082, £ out = 3.654, 3.984, 5.086, for 
ft = 0, ft = 0.5ft c = 0.041, ft = ft c = 0.082, respectively. In Fig.6 we show the 
behavior of M 3 / 2 (po)|a, for different values of ft m = 0, ft m = 0.5/3 c , ft m = ft c , 
for which M3/2 = 2.714, 3.081, 3.622, at po = 1, respectively. 



5 Application to the Local galaxy cluster 

The question about the importance of DE on the dynamics of the Local Cluster 
(LC) was considered by Chernin (2001, 2008) and Chernin et al. (2009). Simple 
estimations confirm this importance and, for the presently accepted values of 
the DE density p v = (0.72 ± 0.03) • 10~ 29 g/cm 3 , the mass of the local group, 
including the dark mater input, is between Mlc ~ 3.5 • 10 12 M Q , according to 
Chernin et al. (2009), and Mlc ~ 1-3 • 10 12 M Q , according to Karachentsev et 
al. (2006). The radius Rlc of the LC is known even worse. It may be estimated 
by measuring the velocity dispersion vt of galaxies in LC and by the application 
of the virial theorem, so that Rlc ~ \J (GMlc Jvt)- The velocity dispersion 
of the galaxies in LC, estimated by Karachentsev et al. (2006), resulted equal 
to v t — 63 km/s, very close to the value of the local Hubble constant H = 68 
km/s/Mpc (Karachentsev et al., 2006). 
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Similarity of this values indicates the great difficulties in dividing the mea- 
sured velocities between regular and chaotic components. In addition, the in- 
dcfinitcness increases because of the unknown level of the anisotropy in the 
velocity distribution. So, the radius of the LC may be estimated as Rlc — 
{GMlc /vf) — (1.5 -7- 4) Mpc, with a very high error box, which we cannot 
estimate properly. Chernin et al. (2009) give different intervals for the mass 
and radius of LC, without using the fixed value of Vt, but considering the whole 
picture of the distribution of the galaxy velocities. They obtained 

1.2 < Mlc < 3.7 • 10 12 Mq and l.K i4 < 1.6 Mpc. (51) 

It is important to note that Chernin et al. (2009) identified the radius Rlc 
of the LC with the radius R v of the zero-gravity force, which is identical with 
the one corresponding to our critical model with /3 — f3 c , in which the average 
matter density is equal to 2p v , as we can see from ([31]) . These estimations, while 
having very big observational errors, indicate the importance of the presently 
accepted value of DE density on the structure and dynamics of the outer parts 
of LC, and its vicinity. 

Polytropic solutions with DE are not well appropriate for describing the 
LC, being the most mass concentrated in two giant objects: our Galaxy and 
M31 (Andromeda). The polytropic model may have better application to the 
rich galactic clusters, where the mass is more uniformly distributed among a 
large number of galaxies and the distribution of matter could be approximated 
by an averaged polytropic distribution. In this aspect, it is necessary to make 
more extended measurements of the galaxy velocities in the clusters and of their 
distribution over the cluster extension. 

6 Conclusions 

The density of DE, measured from SN la distributions and spectra of CMB 
fluctuations, imply the necessity to take into account such a contribution in 
calculations of the structure of galaxy clusters. We developed these calculations 
by considering simple polytropic models, in which it is possible to see how the 
model changes with increasing of the influence of DE. 

Three different values of n have been taken into account. Due to the presence 
of DE, only for central densities larger than a critical value p c , depending of n, 
static solutions are found. We have derived the virial theorem for equilibrium 
configurations in presence of DE, finding expressions between different sorts of 
the energy in polytropic models. 

We have analyzed the stability of equilibrium configurations in presence of 
DE using energetic method. It is shown, that DE increases the instability of 
the equilibrium configurations, working in the same direction as the influence of 
the general relativistic corrections. The presence of DE decrease the the level 
of relativity of a warm dark matter (WDM), at which it is still possible to form 
gravitationally bound configurations from WDM. 
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Finally, the observational indcfinitcncss in the parameters of LC does not 
permit to draw definite conclusions about the level of DE influence but, without 
any doubt, it indicates the dynamic importance of DE in the scale of the galaxy 
clusters. 
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-0.2- P = 

Figure 1: The density distribution for configurations at n = 1 with j3 — 0, /3 = 
0.5/3 c , /3 = /3 C . The curves are marked with the values of f3. The nonphysical 
solution at /3 = 1.5/3 C , which has not an outer boundary, is given by the dash- 
dot line. The nonphysical parts of the solutions at (3 < f3 c , behind the outer 
boundary, are given by the dash lines. The solutions, at large £, asymptotically 
approach the horizontal line 9\ = (3. 
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Figure 2: Nondimensional mass Mi of the equilibrium polytropic configurations 
at n = 1 as a function of the nondimensional central density p n , for different 
values of f3i n - The cosmological constant A is the same along each curve. The 
curves at /3, n ^ are limited by the configuration with (3 — (3 C . 



16 




Figure 3: Same as in Fig.l, for n = 3. The solutions asymptotically approach, 
at large £, the horizontal line 6*3 = /3 1 / 3 . 
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Figure 4: Same as in Fig. 2, for n = 3. 
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Figure 5: Same as in Fig.l, for n = 3/2. The solutions asymptotically approach, 
at large £, the horizontal line 9 3 / 2 = /? 2 / 3 . 
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